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Abstract 

We study in this paper the integration of elastic and streaming traffic on a same 
link in an IP network. We are specifically interested in the computation of the mean 
bit rate obtained by a data transfer. For this purpose, we consider that the bit rate 
offered by streaming traffic is low, of the order of magnitude of a small parameter 
s <C 1 and related to an auxiliary stationary Markovian process (X(t)). Under the 
assumption that data transfers are exponentially distributed, arrive according to a 
Poisson process, and share the available bandwidth according to the ideal processor 
sharing discipline, we derive the mean bit rate of a data transfer as a power series 
expansion in e. Since the system can be described by means of an M/M/l queue with 
a time- varying server rate, which depends upon the parameter e and process (X(t)), 
the key issue is to compute an expansion of the area swept under the occupation 
process of this queue in a busy period. We obtain closed formulas for the power 
series expansion in e of the mean bit rate, which allow us to verify the validity of 
the so-called reduced service rate at the first order. The second order term yields 
more insight into the negative impact of the variability of streaming flows. 

Key words: variable M/M/l queue, perturbation theory, processor-sharing 



1 Introduction 



The emergence of the Internet as the universal multi-service network raises 
major traffic engineering problems, in particular with regard to the coexistence 
on the same transmission links of real time and data services. As a matter 
of fact, these two types of services have different requirements in terms of 
transfer delay and loss, data transmission being very sensitive to packet loss 
but relatively tolerant to delay whereas real time services have strict transfer 
delay constraints. While classical data transfers are usually controlled by TCP 
(Transmission Control Protocol), which aims at achieving a fair bandwidth 
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allocation at a bottleneck link (see Massoulie and Roberts (Jl|) for a discussion 
on modeling TCP at the flow level and processor sharing), real time services 
most of the time are supported by the unreliable UDP protocol, even if some 
transmission control can be performed by upper layers (e.g., RTCP). Real 
time services thus reduce the transmission capacity for data transfers. 

This problem has been addressed by Delcoigne et al. (@), where stochastic 
bounds have been obtained for the bit rate seen by a TCP data transfer, 
when elastic traffic and unresponsive streaming flows are multiplexed on a 
same link (see also Bonald and Proutiere (0)). From a theoretical point of 
view, this problem can be seen as the analysis of a priority system, where 
streaming flows have priority over data traffic. In this context, a usual ap- 
proximation (referred to as Reduced Service Rate, RSR) consists of assuming 
that everything happens as if the service rate for data were reduced by the 
mean bit rate offered by streaming flows. This approximation has been inves- 
tigated for the number of active data flows by Antunes et al. (jH), when the 
load offered by streaming flows is very small. We note that the same kind of 
problem has been addressed in the technical literature by Nunez-Queija and 
Boxma (jD) in the context of ABR service in ATM networks and more recently 
by Nunez-Queija (0; 0) via matrix analysis for systems described by means 
of quasi birth and death processes. In a similar context, Nunez-Queija et al. 
(0) use a perturbation technique for studying a priority system, where prior- 
ity traffic offers a small load. Although the systems considered in this paper 
are quite general, no explicit expressions for the terms of the expansions are 
provided. Finally, note that systems with different speeds are also of interest 
for analyzing the coexistence of different traffic types (J9|). 

In this paper, we investigate the mean bit rate obtained by a data transfer 
when elastic traffic and unresponsive streaming flows are multiplexed on a 
same transmission link. Along the same line of investigations as Antunes et 
al. (jij), because of the real difficulty of the problems, the mean load offered by 
streaming flows is supposed to be very small (controlled by a parameter e < 1) 
and a perturbation analysis for the analysis of the mean bit rate is done. It is 
assumed that elastic flows arrive according to a Poisson process and share the 
available bandwidth according to the processor sharing discipline. In addition, 
to simplify the computations, we assume that the service time required by 
data transfers is exponentially distributed with parameter \x. Thus, we have 
to deal with an M/M/l queue with a time- varying server rate, which depends 
upon the instantaneous number of active streaming flows. The exponential 
distribution of the service time is of course not realistic in practice, but as it 
will be seen later, the mathematical analysis is already very difficult in the 
case so that our work should be considered as a first step in this domain. 

To compute the mean bit rate of a data transfer, we consider the quantity 
A = Jq L(s)ds, where B is the length of the busy period and L(t) is the 
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number of customers at time t in the M/M/l queue under consideration. The 
quantity A is equal to the cumulative waiting time in the M/M/l queue and 
also represents the amount of data served during a busy period. In the case of 
the M/M/l PS sharing queue, if E(cf) represents the mean bit rate obtained 
by a data transfer, we have E(A) = E(N)E(S) /E(d) = E(B)/E(d), where 
E(S) is the mean service time (equal to 1///) and E(N) is the mean number 
of customers served in a busy period. Thus, the computation of E(A) allows 
us to estimate E(d), since the quantity E(B) has been computed by Antunes 
et al. (4) as a power series expansion of e. Note that in the case of a classical 
M/M/l queue, we have E(B) = l/(/i(l-p)) and E(A) = l/(/i(l-p) 2 ), which 
yields E(d) = (1 — p), where p is the offered load. This is the classical result 
for an M/G/l PS queue, which states that the mean bit rate obtained by a 
data transfer is (1 — p) times the server rate (taken as unity in this paper); 
see Massoulie and Roberts (|l|). 

In this paper, we derive a power series expansion in e of the quantity E(A) in 
the case of an M/M/l queue, whose server rate is modulated by an auxiliary 
process (X(t)). We specifically assume that the server rate at time t is p + 
ep(X(t)) for some function p satisfying regularity assumptions, the process 
(X(t)) being stationary, ergodic, and Markovian. The objective of this paper 
is, first to check the validity of the RSR approximation, which claims that 
everything happens as if the server rate were frozen at the value //+eE[p(X(0))] 
and, second, to get some qualitative insight on the impact of the variability 
of streaming flows on elastic traffic. 

The organization of this paper is as follows: The model is described in Sec- 
tion 2, where the main result concerns the power series expansion in e of E(^4), 
which is the key quantity for computing the mean bit rate of a data transfer. 
In Section 3, the main result is applied to obtain the expansion of the mean 
bit rate of a data transfer and to analyze different special cases. Some con- 
cluding remarks are presented in Section 4. The quite technical proof of the 
main result is sketched in the Appendix. 



2 Model description 

Throughout this paper we consider a stable M/M/l queue with arrival rate 
A and service rate p; the load p = \/p < 1. Let L(t) denote the number of 
customers at time t. The invariant distribution ir of (L(t)) is geometrically 
distributed with parameter p. 

Let B denote the duration of a busy period starting with one customer, that 
is, B = inf{s > : L(s) = 0}, given L(0) = 1. For x > 1, let B x denote the 
duration of a busy period starting with x customers. Note that £>i d =' B. In 
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the following, when the variables B, B\ and B[ are used in the same expression, 
they are assumed to be independent with the same distribution as B. 

The quantity A defined in the Introduction represents the area swept under the 
occupation process in a busy period. When several busy cycles are considered, 
the notation Ab will be used to indicate that the area is calculated for the 
corresponding busy period of length B. By definition, the relation A > B 
holds, the excess will be denoted by A = Ab = Ab — B. This queue will be 
referred to as the standard queue denoted, for short, by S-Queue. 

Streaming flows impact data transfers by reducing the amount of available 
bandwidth. This situation is described by introducing an M/M/l queue with 
arrival rate A and varying service rate driven by an ergodic Markov process 
(X(t)) taking values in a state space S. Typically, the state space of the envi- 
ronment is a finite, countable set when (X(t)) is a Markov Modulated Poisson 
Process or S = R in the case of a diffusion, for instance an Ornstein-Uhlenbeck 
process (see Fricker et al. |lfll )). The invariant measure of the process (X(t)) 
is denoted by v. The Markovian notation ~K X (-) will refer only to the initial 
state x of the Markov process (X(t)). 

Let L £ (t) be the number of customers of the queue at time t. The process 
(L e (t), X(t)) is a Markov process. If X(t) = x and L(t) = n > 0, then the 
service rate is given by p + ep{x) for some function p(x) on the state space 
of the environment S and some small parameter e > 0. For t > 0, let us 
define the quantities p + (t) = max(p(t),0) and p~(t) = max(— p(t),0) so that 
pit) = p + (t) —p~(t). At time t, the additional capacity is therefore ep + (X(t)) 
and ep~(X(t)) is the capacity lost. 

In the rest of this paper, we make the two following assumptions: 
the function p(x) is bounded (Hi) and A + e sup(|p(x)| : x G S) < p, (H2). 
If assumption (H 2 ) holds (it is a suflifient but not necessary condition) then 
the queue is stable and, in this case, the duration B £ of a busy period starting 
with one customer, B e = inf{s > : L £ (s) = | £ e (0) = 1}, is a.s. finite. The 
queue with time- varying service rate as defined above will be referred to as the 
perturbated queue, denoted, for short, by P-Queue. The case e = obviously 
corresponds to the S-Queue. The area for the perturbated queue over a busy 
cycle is defined as 

A £ = / L e (s)ds. 
Jo 

The basic idea of the perturbation analysis carried out in this paper for the 
quantity A £ defined by the above equation is to construct a coupling between 
the busy periods of the processes (L £ (t)) and (L(t)). Provided that for both 
queues the arrival process is the same Poisson process with parameter A, we 
add and remove departures as follows. 
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Additional departures. When p + (X(t)) > 0, there is additional capacity 
when compared with the S-Queue and more departures can take place. 
These additional departures are counted by means of a point Process Af + = 
(tf), with < tf < tf < ■ ■ ■ , which is a non-homogeneous Poisson process 
on M + with intensity given by t —>■ ep + (X(t)). Conditionally on (X(t)), 
the number of points of jV + in the interval [a, b] is Poisson with parameter 
£ JaP + (X(s)) ds. In particular the distribution of the location tf > of the 
first point of Af + after is given, for x > 0, by 



P(t+ >x)= P(7V+([0,x]) = 0) = E 



exp 



p + (X(s))ds 



Removing Departures On the other hand, when p~(X(t)) > 0, the server 
rate is smaller than in the S-queue. Let jV^ = (U), a Poisson process with 
intensity \x on R + which represents the non- decreasing sequence of instants 
when the customer of the S-queue (if not empty) may leave the queue. We 
denote by M~ the point process obtained as follows: For s > 0, a point at s 
of the Poisson process j\f^ is a point of Af~ with probability ep~(X(s))/ '//. 
(Note that this number is < 1 by assumption (H2) ) The point process j\f~ 
is Poisson with intensity s — > ep~(X(s)). A point of j\f~ is called a marked 
departure. The points of J\f~ are denoted by < t{ < t^> < • • • < t~ < ■ ■ ■ . 
By definition, 



P(t- > x) = E I J| 1 - 

\ti,ti<X \ 



x>0. 



The processes defined above are non homogeneous Poisson processes (see 
Grandell ffllh for an account on this topic). The main result of this paper 
is the expansion in power series of e up to the second order of E( v 4 £ ). 



Theorem 1 The second order expansion of the area swept under the occupa- 
tion process of the perturbated queue during a busy period is given by 



e a<) — e (a) - £ m^m^±A _ , 2(a+ + a _ + a±) + 0(fi2)> 

v ' (m - A) 3 



(3) 



where the coefficients a + , a_, a± are defined below by Equations (A. 4) , (A. 5), 
and (A. 6), respectively. 



3 Applications 



In this section, as an application of Theorem 1, we evaluate the mean bit rate 
E(d e ) obtained by an elastic data transfer when the server rate is perturbated 
by the presence of streaming flows (coming through the term ep(X(t)) in the 
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service rate of the perturbated M/M/l queue, equal to \i + ep(X(t))). As 
mentioned in the Introduction, we have E(^4. £ ) = E(B £ ) /E(d e ) . The average of 
the duration B £ of the corresponding busy period has been studied by Antunes 
et al. (4) and can be expanded in power series of e as follows. 

Theorem 2 The expansion ofE(B £ ), the mean duration of a busy period, is 
given by E(B £ ) = 1/0 - A) - eE u (p(X (0))) / (p - A) 2 + (6_ - b+)e 2 + o(e 2 ). 
where b + and 6_ are given by, with the notation of Proposition 8, 

b + = ~E ( J\b - v) E„ (p + (X(0))p + (X(v))) dv^j 



E EE/ P + (X(u))p-(X(Dl))du\, (4) 



p\l-p) 



A=l j=l 



l ( ( N r B+Bi 



p\l-p) 

1 



AT AT' 

+; e IEEp1 i ( I) '))p"W 5 + 4)) 1 1 • (5) 



By using Equations (4) and (5) and Theorem 1, straightforward computations 
show that the quantity E(d £ ) can then expanded in power series of e as follows. 



Proposition 1 The mean bit rate of an elastic data transfer can be expanded 
in power series of e as 

m)=l _ p+ e^Mmi E+ ^ +0( ^ (6) 

n 

where the coefficient c is given by 

c = E.(p(X(0))) 2 ^ (1 l ^ ) + p(l - pf ((1 - p){a + + a_ + a±) + 6- - 6+) , 

^ ^ (7) 
the quantities a + , a_ ; a±, b + and 6_ ieinj defined by equations (A. 4), (A. 5), 
(A. 6), (4) and (5), respectively. 

From Equation (6), we immediately deduce that as far as the first order term is 
concerned, the RSR approximation is valid: for an M/M/l queue with service 
rate p + eE u (p(X(0))), the mean bit rate denoted by d obtained by a customer 
is given by E(d) = 1 - A/0 + eE u \p(X(0))]) = 1 - p + £pE u [p(X(0))]/ p + o(e) 
Unfortunately, the coefficient c defined by Equation (7) intricately depends 
upon the correlation structure of the modulating process (X(t)) and the dy- 
namics of the M/M/l queue. Because of this complexity, three cases of prac- 
tical interest in the following are considered: non-positive perturbation func- 
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tions, non-negative perturbation functions, and special environments (namely, 
fast and slow environments). 

3.1 Non-positive Perturbation Functions 

We assume in this section that the perturbation function is non-positive so 
that the environment uses a part of the capacity of the M/M/l queue with 
constant service rate \i. This application is motivated by the following practical 
situation. Coming back to the coexistence of elastic and streaming traffic in 
the Internet, assume that priority is given to streaming traffic in a buffer of 
a router. The bandwidth available for non-priority traffic is the transmission 
link reduced by the rate of streaming traffic. Denoting by er(X(t)) the rate 
of streaming traffic at time t (for instance e may represent the peak rate 
of a streaming flow and r(X(t)) the number of such flows active at time 
t), the service rate available for non-priority traffic is /i — er(X(t)). Setting 
p(x) = —r(x), the function p(x) is non-positive. 

Proposition 2 When p + = 0, if d is the mean bit rate of the M/M/l queue 
with service rate fj, + eE[p(X(0))] ; then with the notation of Proposition 8, 

hmiE(4-rf) = -^^E( £ C P (XP,-A))] 

fc fJ> \ 1 Si^AS AT I 



n _ n \ 3 ( N N ' \ 

V—y-E ee^-a+d;)(bi-d;) , (8) 

^ \i=l 3=1 ) 



where the function C p (u) is defined for u > by 

C p (u) = E„ [p(X(0))p(X(u))} - E v [p(X(0))} 2 (9) 

and is, up to the factor e 2 , the auto- covariance function of the variable capacity 
of the perturbated queue. 

The above result shows that if the process (X(t)) is positively correlated, i.e. 
Cp(') > 0, then lim e ^oE (d £ — dtj /e 2 < 0. The environment has therefore a 
negative impact on the performances of the system in this case. 

Proof. The terms a+, a± and b + in the coefficient c defined by Equation (7) 
are equal to 0. In addition, one easily checks that 



(i _ p)a _ + b _ = _i_ E ( £ p -(x(A)K(x(/J J )) 

K l<i<j<N 
N N' 



VP 



p Vi=ifc=i 



One concludes by using Expansion (6). 
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3.2 Non-negative Perturbation Functions 

It is assumed in this section that p~ = 0. We have the following result, which is 
the analogue of Proposition 2 for this case. Contrary to the above proposition, 
the expansion has a more explicit expression. Its (straightforward) proof is 
omitted. 

Proposition 3 When p~ = ; with the same notations as in Proposition 2. 



An integration by parts and some calculations give the following corollary. 

Corollary 3 When the correlation function of the environment is exponen- 
tially decreasing, i.e., when C p (x) = Var[p(X(0))] e~ ax for all x > and some 
a > 0, then 



with the convention that, if Z is some integrable non-negative random variable, 
the density on M.+ of the variable Z* is defined as 

¥(Z* >x) =F(Z >x)/E(Z), 

for x>0, and Z** stands for (Z*)*, and Z*** for (Z**)*. 

When a is small, the right hand side of Equation (11) is equivalent to the 
quantity — 2p//x 2 < 0. This shows that non-negative perturbation functions 
have a negative impact at the second order on the mean bit rate of elastic 
data transfers. 

Remark. This result could maybe be extended to a perturbated function 
with a non-constant sign but it is out of reach because of the complexity of 
the secund order term. 

3.3 Fast and slow Environments 

The performance of the system in two limit regimes, called fast and slow envi- 
ronments, are now evaluated. These regimes are very useful, since performance 
in the limit regimes is insensitive and only depends on appropriately defined 
parameters. Such a situation has also been analyzed by Delcoigne et al. (@) 
through stochastic bounds. 




(10) 
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The environment is scaled by a factor a > 0, such that at time t the environ- 
ment is supposed to be X(at). The behavior when a goes to infinity and zero 
is investigated. 

When the parameter a is very large, the environment process approximately 
averages the capacity of the variable queue. For a large a and for t and h > 0, 
the total service capacity available during t and t + h is given by 

fxh + e f t+h p(X(au)) du d = ph + e- p{X{u)) du~(p,+ eE\p(X (0))})h 
Jt a Jo 

using the stationarity of {X(t)) and the ergodic theorem. Thus, when a tends 
to infinity, the variations completely vanish and the service rate reduces to a 
constant. 

On the other hand, for small values of a, the environment process remains 
almost constant over the busy period of the P-Queue. As a goes to 0, the 
variation disappears and the environment is frozen in the initial state of the 
process: the service rate is constant and equal to /i + ep(X(0)). 

This intuitive picture is rigorously established in the next proposition. In the 
following a general perturbation function p is considered together with some 
stationary Markov process {X(t)) with invariant probability distribution v. It 
is assumed that it verifies a mixing condition such as 

lim \E[f{X(0)g(X(t))]-E v {f)E„{g)\=0 (12) 

t — >+oo 

for any Borelian bounded functions / and g on the state space S. Note that this 
condition is not restrictive in general since it is true for any ergodic Markov 
process with a countable (or finite) state space or for any diffusion on M. d . 

Under the above assumptions, we have the following result; the proof relies 
on the use of the mixing condition (12) and can be found in the paper by 
Antunes et al. |l2l). 

Proposition 4 When the environment is given by (X(at)) and Relation (12) 
holds, then when e tends to 0, E(d e ) — E(<f) = ^/(a)e 2 + o(e 2 ) where 

.. T , . E,(p(X(0))) 2 JV T . . E u (p(X(0)) 2 ) 
km Via) = -p — n and limf a = -p „ v " ' . 

The fast and slow environment provide an explicit estimate of the second 
order term, where the slow environment yields a worst performance of the 
perturbated queue. It is not clear that these limit regimes give a lower and 
upper bound of the performance of the queue. 
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4 Conclusion 



We have investigated in this paper the impact on the performance of elastic 
data transfers of the presence of streaming flows, when both kinds of traffic 
are multiplexed on a same link of an IP network. By assuming that the pertur- 
bation due to streaming flows is of small magnitude, a perturbation analysis 
can be performed in order to obtain explicit results for the mean bit rate 
achieved by a data transfer, under the assumption that elastic streams share 
the available bandwidth according to the processor sharing discipline. It turns 
out that at the first order, the so-called RSR approximation is valid. This is 
not the case for the second order, for which the variability of streaming flows 
seem to have a negative impact, at least for the three cases examined here. 

Further investigations are needed in order to estimate the degradation suffered 
by data transfers at the second order. The perturbation analysis carried out 
in this paper is possible, because we have assumed that streaming flows offer 
a very small contribution to the total load. When this is not the case, new 
tools have to be developed to estimate the quality of data transfers. 



A Appendix: Proof of Theorem 1 



Since the derivation of the expansion of the average area is somewhat technical, 
we begin with the simplest case, which is the first order expansion. In the 
following, we set A £ = E(A) - E(A £ ). 



A.l First order term 



We first consider an additional departure. Define £ + = {tf < B,t 1 > tf + 
B L / t +\_ 1 \, where B L ( t +^_ 1 is the duration between the tf and the first time 
when the S-Queue has one customer (see Figure A.l-(b)). On this event, an 
additional departure is added and the busy period of the P-Queue finishes 
before a departure is canceled. For the first order term of the expansion of the 
mean value of A £ on the event S + , we only need to consider the case where 
there is only one additional departure during the busy period of the P-Queue. 
The probability that two additional jumps occur in the same busy period is 
of the order of magnitude of e 2 since the intensity of the associated Poisson 
process is proportional to e. 
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Lemma 5 In the case of a single additional departure 



E(A .- l£+) ^ww±^ + 0(E) . 

Proof. The difference A e on the event {tf < B, t£ > tf + > 
tf + #£,( t +)_i} is the sum of two disjoint areas (see Figure A.l-(a)). The first 
one is given by the distance between tf and the end of the busy period of 
the S-Queue. By the strong Markov property at the stopping time B £ , con- 



B 



L(t 



ditionally on the event {tf < B,tf > tf + B L ^ t +- ) _ 1 ,t 1 > t\ 

the S-Queue starts at time B £ an independent busy period with one customer 
(with duration Bi). The second area of A £ is then given by the area of the 
sub-busy periods in B ± , i.e. periods when £— is > 1 in the second b.p. 



P-Queue 
S-Queue 




Fig. A.l. (a) One additional departure — (b) One marked departure 



It follows that 



E (A £ l, + ) = E {*h^^ + B^_^ + B^ + o(e) = 

+ E(A Bl - B 1 )F(tf <B,tf>tf + B L{tt) _ x ,tl >tf + B L{tt) _f) 

= F(tf < B)E(A Bl — Bi) + E ((B - tf)l {t +< B} ) + o(e) 



Equation (1) and the boundedness of p give that 



P(t+ < S) = 1 - E 



exp — e 



l B p + {X{s)) 
Jo 



ds 



eE 



f B p + (X(s)) 
Jo 



ds 



+ o(e) = £E(5)E p + (X(0)) 



+ £ 



-E 



p + (X(0)) +o(e) 



by independence between B and (X(t)) and by the stationarity of (X(t)). 



11 



Similarly, 



e((S-^)% b} )=e|| ep + (X(u))e- £ i:^ x ^ ds (B-u)du) 

= eE[p+(X(0))]^fl + o(e) 

Since E (B 2 ) = 2/{i 2 (l - p) 3 and E(A Bl ) = ^/{^ - A) 2 (see for instance 
standard books such as Cohen (flif))). The lemma is proved. □ 

We now turn to the case, when there is one removed departure. On the event 
£_ = {ti < B, B + B 1 < tf}, a marked departure occurs and no departures 
are added before the completion of the busy period B\. We derive the first 
order expansion of the mean value of A £ on 

Assume that there is only one marked departure and no additional jumps 
during the busy period of the P-Queue. In this case, at the end of the busy 
period of the S-Queue, the P-Queue has one customer and the difference is 
the distance between t\ and the end of the busy period B. At time B, the 
P-Queue starts a busy period with one customer and provided that there are 
no marked and additional departures during (B,B + B e ), the difference has 
the same distribution as the area of a busy period B\ of the standard queue 
(see Figure A.l-(b)). 

Lemma 6 In the case of a single marked departure 

E(4%) .. E !tfflM ++) . 



Proof. By using the same arguments as before, one obtains the relation 

E(A S 1 £ _) =E(Atl {t -^ t -> B+Bi4 > B+Bi} ) +o(e) 

= -E ((B -f{ + A Bl )l {tI < B , t -> B+Bl)tt > B+Bl} ) + o(e) (A.l) 

Hence, E(A £ 1^_) = -E ((5 - + A Bl )t {tI < Bit -> B+Bljt +> B+Bl} ) + o(e), so 
that E(A £ 1 £ _) = -E((B - t^)l {t -< B} ) - E{A Bl )F(ti < B) + o(e). 

To estimate F(ti < B), let (Di) denote the sequence of departures times and 
the number of customers served during the busy period of length B, then 
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Equation (2) gives the identity 

P(f - < B) = E (£*p-W» n f i - 

\i=l ^ j=l V Z 1 / 

= £ E (l>-(*(A))l + o{e) = -E(N)M\p-(X(D 1 ))} + o(e) 

by stationarity of (X(t)) and Wald's Formula with E(iV) = 1/ (1—p). Similarly, 

K((B-tr)l„ r < B) ) = iE fep"WA))(B-A))+o(e) 

= £ E 1 p- , (l-W))] 

where _D = I^ili -A is the sum of the departures in the busy period of the S- 
Queue. Using the fact that E(D) = /j 2 / (p — A) 3 , E(NB) = (1 +p)/(>(l - p) 3 ) 
and E(A) = p/(p — A) 2 , the result is proved. □ 

Combining Lemmas 5 and 6 yields the first order term indicated in Theorem 1. 
A. 2 Second order term 



To compute the second order term in the power series expansion in e of E (A 
three cases have to be considered: 

- tf<B £ or t2<B £ : one or two additional departures occur in a busy period; 
— t\ < B £ or t% < B £ : one or two departures are canceled; 

- ti < B £ and < B £ : one additional departure takes place and another 
one is canceled. 

It is not difficult to show that any event involving a third jump yields a term of 
the order e 3 in the expansion of the mean bit rate. Due to the space constraints, 
a part of the expansion is proved. The complete proofs of the expansion can 
be found in Antunes et al. (II 



In the case that two additional departures occur during B £ , the difference 
between the areas of the busy periods due to the first additional jump is given 
by the two first terms on the right hand side of the following equation 

E ( A£ V<^2 + <^+S Mt + ) _ 1 }J = E [( B ~*^ 1 {^<B,4<^+-B i(t + ) _ 1 }) 

+ E(^ Bl )P(t+ < B,t+ < ^ + 5 i(t + ) _ 1 )+Efs L(t + ) _ 1 l {t + 

<B,t%<t++B + _ } J 

+ E(,4bi )P(t+ <B,t+<tt + B L{t{yi ) + o(e 2 ), (A.2) 
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which follows by the same arguments stated for only one additional departure. 



P-Queue 

S-Queue 




a B 



Fig. A. 2. Two additional departures 

Due to the second additional jump, the difference A £ increases by the sum 
of two disjoint areas. The first one is given by B L , t +-,_ 1 and represents the 
distance from the second additional jump until the first time the S-Queue 
with less than one customer is empty. Note that by conditioning on the event 
{tt < B,tt <tt + ^L(t+)_i}> & t the stopping time B £ , a new busy period B[ 
starts with the same distribution as B. Thus, the second area is given by the 
area of sub-busy periods in \wB\ (periods where t — > L(t) is > 1) in B[, hence 

E ( A ^+) = E ( A£1 ^<^>^ +BMt , ) _ 1 ,r>^H- Bi(t , ) _ l} ) 

+ E { A %ti<B4<4 + B L(tt) _j)+o(e 2 ) 

=E ((B - t+)l { ,+< B} )+E (^ (t2+) -il { ^<s4<^+^ (t + ) _ 1 }) +HA Bl )ntt<B) 
+E(A B[ )¥(tt <B,t+<tt + #L( t +)-i) + K(e) + o(s 2 ). (A.3) 

where K(e) is a term which is not expressed here for sake of simplicity. The fol- 
lowing proposition gives the expansion of some of the terms of Equation (A.3). 
The other expansions are done much in the same way (with various compli- 
cations). 

Proposition 7 The following expansions hold 



F(tt<B)= 

Eb+(X(0))] 



// — A 



-e 2 E f J (B-v) E [p + (X(0))p + (X(v)) 



dv +o(e 2 ), 



ntt<BA<tt+B Lit+i) _ l ) 

=e 2 pE 



J\b-v)E 



p + (X(0))p + (X(v))] dv^j +o(e 2 ) 
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Proof. Since ¥(tf < B) = E (l - cxp (-e j* B p + {X(s)) dsj), the expansion in 
power series of e has the first term E(p + (X (0))) / (fi — A) and second term 



y|| B j) + (I(s))(ls| )=E^J <B E(p + (X(0))p + (X(v-u)))dudv 



B 







B-v)E \p + (X(0))p + (X(v)) 



dv 



by stationarity of the process (X(t)). The first expansion is proved. 

The event {ti < B ,t\ < tf + B L ^_- i } occurs only when tf and t\ are in 
a sub-busy period [sj_i + sj, for some i G {1, . . . , if}, (a period where 
L is always > 1). The variables E^ are i.i.d. exponential with parameter A, 
B\ — Si — Sj_i — x has the same distribution as B and H is geometri- 
cally distributed with parameter A/ (A + //). The probability that the first two 
additional jumps are in the i-th sub-busy period, is 

E {J* +e ep + {X{u))e-^> + ^ ds (l - e -J> + W*))<^ ^ 

= e 2 e( T p + (X(u)) r p + (X(s))dsdu) +o(e 2 ) 

= 5 2 E f / E„ (p+(X(0))p + (X(t; - «))) du) + o(e 2 ), 

\J0<u<v<B v ' J 

which gives the second expansion. □ 
The complete expansion is now detailed. 

Proposition 8 The coefficients of e 2 in the expansion of E(A e ) are given by 

B r > 

dv 



a + = -—j-^—^E \J o (B-v)E [p+(X(0))p + (X(v)) 

1 ~ P -E [J\b - u) 2 E [p + (X(0))p + (X(u))] du) , (A.4) 



-E £ p-{X{D,^)p-{X{D 3 )) 



/A 1 - P) \l<z^<N 

^ (e&IWWB + m (b, - D>. + j , (A.5) 
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« ± = ±* [Y.l B r(x(D, ))p %x { s ))d s 

+ > + ,)) ( Bl - . + ^) *) 

- ^ (e E /* («))p"(* (A* )) ((^Ta)5 + ^ - ") <M < A ' 6 > 

where H is geometric distributed with parameter A/(/i + A), (iVj, .D^, . . . -Dyv ) 
denotes respectively the number of departures and the departures times in a 
busy period B\, and A4 — B\ + E + J2k=i+i(^k + -Bi) where (Ei) are i.i.d 
is exponentially distributed with parameter /i + A and (B\) are i.i.d with the 
same distribution as B. 
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